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ABSTRACT 

The highest velocity neutron stars establish stringent constraints on natal kicks, asymmetries in 
supernova core collapse, and the evolution of close binary systems. Here we present the first results of 
a long-term pulsar astrometry program using the VLBA. We measure a proper motion and parallax 
for the pulsar B1508+55, leading to model-independent estimates of its distance (2.37+q'2q kpc) and 



transverse velocity (lOSS+gQ 3 km s 



J^gQ 3 km s _1 ), the highest velocity directly measured for a neutron star. We 
trace the pulsar back from its present Galactic latitude of 52.3° to a birth site in the Galactic plane 
near the Cyg OB associations, and find that it will inevitably escape the Galaxy. Binary disruption 
alone is insufficient to impart the required birth velocity, and a natal kick is indicated. A composite 
scenario including a large kick along with binary disruption can plausibly account for the high velocity. 
Subject headings: astrometry — stars: neutron — pulsars: individual (B1508+55) — stars: kinematics 
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1. INTRODUCTION 

The high velocities of young radio pulsars relative to 
their progenitor stellar populations have been readily 
apparent, both in their proper motions, and through 
the bow shock nebulae produced in some cases. Mod- 
els to explain these high neutron star (NS) velocities 
have invoked the disrupt ion o f close binary system s 
llBlaauwlll96it IGott et al]ll97(l llben fc Tutiika\Tll99fi|) . 
natal kicks from asymmetric supernova (SN) explosions 
(Shkl ovskiilll970tlDewev fc Cordell987l) . and post-natal 
acceleration through the ele ctromagnetic rocket effect 
ijHarrison h Tademarul Il975j) . As the pulsar popula- 
tion velocity distribution has been elucidat ed (see, e.g., 
lArzoumanian et aT]l2002HHobbs et al.ll2005j) . binary dis- 
ruption and the electromagnetic rocket effect have re- 
quired more extreme assumptions to remain viable, and 
natal kicks have emerged as the most plausible mech- 
anism for NSs to gain high velocities. In turn, ob- 
served pulsar velocities have placed stringent constraints 
on SN core collapse mechanisms. In simulations, a na- 
tal kick of the correct order of magnitude has been 
obta ined from hydrodynamic and convective inst abili- 
ties l)Burrows fc HavesTl996tlJanka fc Muellerll99fiD . and 
more exotic mechanisms involving asymmetric neutrino 
emission in the presence of strong magnetic fields (B > 
10 15 G) have also been suggested (jArras fc Laill99fl) . 

At this point, two caveats are worth careful attention. 
First, the current pulsar velocity distribution suffers from 
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severe selection effects. The stellar progenitors of NSs in- 
habit the Galactic plane, as do we, and thus higher veloc- 
ity pulsars spend less time within our detection volume 
than do lower velocity ones. The effect is exacerbated 
by pulsar searches that focus more heavily on the target- 
rich environment of the Galactic plane. The extremely 
successful Parkes Multibeam survey, for e xample, spans 
only \b\ < 5° (e.g ^Manchester et all2 001). Thus the ob- 
served velocity distribution is likely to be biased toward 
lower velocities. 

Second, most pulsar distances (and hence veloci- 
ties) are inferred from their pulsar dispersion mea- 
sures (DM = Jq n e dl) and models for the Galac- 
tic distributio n of electron density (n e ) along the 
line of sight llTavlor fc Cordesl ILoM hereafter TC93; 
iCordes fc Lazicl2005l hereafter NE2001) which have sub- 
stantial uncertainties. For example , the Guitar Neb- 
ula pulsar B2224+65 l)Cordes et all 119931) is a show- 
case example of ultra-high transverse velocities, V± « 
1640 (D/l. 9 kpc) km s _1 , but the distance uncertainty 
in the NE2001 model implies that V± < 1000 km s" 1 can- 
not be ruled out. The highest extant model-independent 
estimate of a NS velocity is 630 ± 40 km s~ 4 for 
PSR B1133+16 JBrisken et al.ll200l) . 

Here we present the first hyper- fast pulsar (V > 
1000 km s -1 ) with a model independent distance and 
transverse velocity measurement. PSR B1508-I-55 was 
part of a recently concluded astrometry program with 
the NRAO Very Long Baseline Array (VLBA), in which 
27 pulsars were observed for eight epochs each, span- 
ning 2 years. This Letter serves to introduce the project, 
while comprehensive papers on our astrometric methods 
and parallax results for ^20 pulsars are in preparation 
(W. F. Brisken et al., S. Chatterjee et al.). Below we 
report astrometric results for the distance, velocity and 
birth site of B1508+55, and discuss the implications of 
the high observed velocity. 

2. OBSERVATIONS AND DATA ANALYSIS 

B1508+55 was observed with the VLBA over eight 
epochs spaced by approximately 3 months between 2002 
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Fig. 1. — Left: The motion of PSR B1508+55 in Right Ascension and Declination, with the best fit proper motion and parallax model 
overplotted. The error estimates for each data point are smaller than the size of the points. Right: The parallax signature of PSR B1508+55 
in Right Ascension and Declination, after subtracting the best-fit proper motion from the astrometric positions. Curves corresponding to 
the best fit parallax n = 0.415 mas are overplotted. The inner error bars indicate the random position uncertainties, while the outer error 
bars indicate the net uncertainties (random and systematic, added in quadrature). 



TABLE 1 
Parameters for PSR B1508+55 



Parameter 



Value 



Epoch 

cto 
So 

[la (mas yr _1 ) 
/is (mas yr" 1 ) 
7r (mas) 

D (kpc) 

V x (km s" 1 ) 

l,b 

«,Mb (mas yr -1 ) 

Magnetic Field B (G) 
Spindown Age (yr) 
£>DM ( k P c ) 



2003.0 (MJD 52640) 
15 h 09 m 25?6211 
55°31'32^'331 
-73.606 ± 0.044 
-62.622 ±0.088 
0.415 ±0.037 



q o <y ~\~ . 2 3 
-0.20 

1083 j^g 3 

91.325°, 52.287° 

-6.39, 96.06 



1.95 x 10 12 
2.34 x 10 6 

1.9 (TC93), 0.9 (NE2001) 



Note. — All astrometric parameters are in the 
J2000 coordinate system. Absolute positions are 
referenced to J1510+5702, and have an uncertainty 
~1 mas. in, hi,, and Vx are corrected for differential 
Galactic rotation. 



June and 2004 April. As a trade-off between increasing 
pulsar flux at lower frequencies and improved resolution 
as well as reduced ionospheric effects at higher frequen- 
cies, the observations were conducted between 1.4 and 
1.7 GHz, with four frequency bands of 8 MHz each and 
dual polarization in each band. In order to retain visibil- 
ity phase-coherence, observations were phase-referenced 
by nodding between the target and a primary calibrator 
source (J1510+5702) ~ 1.5° away, with a cycle time of 
120 s on target and 90 s on the calibrator. Residual cali- 
bration errors (primarily from the unmodeled ionosphere 
over 1.5° ) are much reduced by employing in-beam ca l- 
ibration ijFomalont et al J 119991 fChatteriee etafl l200lj) . 



which reduces the angular throw and eliminates the need 
for time interpolation. Thus the observations are finally 
referenced to a faint (~ 7.1 mJy) extragalactic source 
(J151148±554156) only 22.6' away from the pulsar. The 
signal-to-noise ratio (S/N) for the pulsar was boosted by 
gating the correlator on at the expected times of arrival 
of pulses, using current pulse timing solutions obtained 
for each epoch from ongoing observations at the Jodrell 
Bank Observatory. 

Data reduction was conducted in AIPS using a cus- 
tomized pipeline (W. F. Brisken et al. 2005, in prepara- 
tion). Briefly, it included amplitude calibration based on 
system temperatures at each antenna, followed by cali- 
bration of the visibility phases, rates, and delays based 
on the observed visibility phases of the primary calibra- 
tor and the in-beam source. Consistency between epochs 
was ensured by constructing models of the primary and 
in-beam calibrators based on all eight epochs of data, 
and using these models to iteratively calibrate each indi- 
vidual epoch. 

The astrometric positions and position uncertainties 
were fit for a reference position, proper motion, and par- 
allax using a linear least squares algorithm. The ran- 
dom position uncertainties were determined by fitting 
Gaussians to the pulsar images at each epoch and fre- 
quency. Systematic errors were estimated by scaling the 
scatter in astrometric positions at each epoch. The scal- 
ing factor of 1.8 was determined such that the final re- 
duced x 2 f° r the fit was 1.0, with 59 degrees of free- 
dom. The fit results are summarized in Table Q an d 
the proper motion and parallax signatures are plotted in 
Figure^ We find a proper motion fi a = —73.606 ±0.044 
mas yr -1 , fxs — —62.622 ± 0.088 mas yr -1 , and a par- 
allax 7r = 0.415 ± 0.037 mas. We verify the robustness 
of the parallax estimate by omitting each epoch in turn 
from the fit and find consistent results (within ~ la). 

3. DISTANCE, VELOCITY AND BIRTH SITE 



The Hyperfast Pulsar B1508+55 
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From the measured parallax and proper motion, we 
calculate the one-dimensional probability distributions 
for the distance (D) and transverse velocity (Vj_). We 
infer the most probable values and the most compact 
68% confidence intervals as D — 2.37*q'2q kpc and V± = 
1083^gQ 3 km s _1 , the first measured model-independent 
NS velocity > 1000 km s" 1 . 

The pulsar is at a large height above the Galactic plane 
(z = Dsinb ea 2.0 kpc), where DM based distance esti- 
mates are sensitive to the model Galactic electron density 
scale height. Based on the pulsar DM of 19.6 pc cm~ 3 , 
the model distance estimate is 1.9*1 2 kpc fr° m TC93, 
and 1.0 ± 0.2 kpc from NE2001. While TC93 is fortu- 
itously close to the parallax distance, NE2001 produces 
a significant underestimate, implying a much lower elec- 
tron density along part of the line of sight. Such a 
reduction may possibly b e due to voids or "chimneys" 
ijNorman fc Ikeuchill989|) associated with the Cygnus su- 
perbubble. The electron density must then decay rapidly 
after ~ 0.8 kpc along the line of sight. 

While B1508+55 is well above the birth 
scale height for t he G alactic pulsar population 
l)Arzoumanian et alJ |2Q02), its proper motion vec- 
tor (m,Hb = —6.39,96.06 mas yr" 1 after correcting for 
differential Galactic rotation) suggests that it was born 
in the Galactic plane. Given the pulsar period P and 
period derivative P, and assuming that the pulsar was 
born with a rapid initial spin and has been braking 
due to a constant torque solely from magnetic dipole 
radiation 8 (i.e., P oc P 2 ~ n , braking index n = 3), the 
spindown age r = P/2P = 2.34 Myr. 

Tracing the pulsar back in the Galactic potential, we 
find that it could have been born 2.34 Myr ago at 
\z\ < 0.2 kpc for a range of modest (unknown) radial 
velocities between and 300 km s , and sample or- 
bits are shown in Figure The pulsar appears to have 
originated in or around the Cyg OB associations, al- 
though we cannot identify the exact birth cluster due 
to the lack of proper motions and accurate distance esti- 
mates for the clusters. Spindown age is known to be an 
imperfect estimate of the pulsar chronological age (e.g . , 
INaravan fc VivekanandllT98lT: ICordes fc Cnernofl1ll998j) . 
If the true age of the pulsar is significantly larger than 
the spindown age, the pulsar birth location would be sev- 
eral kpc beyond the Cyg OB associations in the Galactic 
plane. However, reasonable values of the pulsar radial 
velocity and a range of pulsar ages (or equivalently, a 
range of braking indices or birth periods) produce a birth 
location in the Galactic plane toward the Cyg OB asso- 
ciations. 

We thus establish a self-consistent narrative for 
B1508+55 that ties together its current location, dis- 
tance, and velocity vector with its spindown age and 
birth site. Born in the Cygnus region, B1508+55 is mov- 
ing rapidly away from the Galactic plane and on course 
to irrevocably escape the Galaxy. For v r = 200 km s -1 , 
the implied three-dimensional NS birth velocity Vns ~ 
1120 km s" 1 . 

4. DISCUSSION: IMPLICATIONS OF THE HIGH VELOCITY 
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8 W hile an estimate for P exists for this pulsar IHobbs et alJ 
2004), it is affected by timing noise and does not give a physically 
meaningful value for the braking index. 
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Fig. 2. — Possible orbits for B1508+55, traced back in the Galac- 
tic potential. The solid dot denotes the current pulsar position and 
the thick solid line the path it has followed for an age r = 2.34 Myr 
(the spindown age, with braking index n = 3) and a radial velocity 
v r = 200 km s . Other possible orbits are shown (thinner lines) 
for t = 4.69 Myr (i.e., n = 2 for a small initial birth period) with 
corresponding v r = —300 km s — 1 , and for r = 1.56 Myr (n = 4) 
with v r = 700 km s _1 . Also indicated are the Cygnus super bubble 
(2 kpc away, dashed ellipse) and the Cyg OB ass ociations (solid 
ellipse s) with positions and extents as tabulated bv lUvamker et all 
12001). The starred symbols are the Galactic SN remnants iden- 
tified in this region JCrccn 2004). The solid horizontal line is the 
Galactic plane and the horizontal dashed lines indicate the pulsar 
birth scale height from Arzoumanian et al. (2002) at the distance 
of the Cygnus superbubble. 



The high birth velocity of B1508+55 constrains evo- 
lutionary scenarios for the pulsar. It could have been 
born in an isolated system, in which case a large natal 
kick is required. Alternatively, it could have formed in 
a binary system that was disrupted by a first or sec- 
ond SN. If so, its companion could have been a longer 
lived (and thus less massive) star, or a compact object. 
However, in order to attain the highest possible veloc- 
ity solely from binary disruption, B1508+55 should have 
formed in a first SN and remained bound to a massive 
He star companion in a tigh t orbit that was subseq uently 
disrupted by the second SN ijlben fc Tutuko\ll 996). The 
net velocity of the pulsar then derives from either one or 
two mass-loss events. Given the steeply declining ini- 
tial mass function for stellar progenitors, He cores are 
expected to have masses < 16M Q . The highest veloc- 
ities are obtained for the smallest possible binary ra- 
dius, which corresponds to the He core filling its Roche 
lobe (i.e., a semi-detached binary). Such binary systems 
would be rapidly circularized due to the interaction be- 
tween the NS and the envelope of the massive companion. 
(However, B1508+55 has a canonical field strength, and 
thus cannot have accreted much material.) For a circu- 
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lar orbit, V NS = \/GM c /a(l + M NS /M C ), where M c is 
the companion mass, a is the orbital radius, and Mns 
is the NS mass (assumed to be 1.4 M ). When such a 
system is disrupted by a symmetric SN explosion of the 
massive companion, the fraction of the orbital velocity 
impa rted to the NS depends on the mass ratio Mns /M c 
(e.g.. lBoersmall96ll). Using estimates for the companion 
star radius iPaczvnskilll971|) and the size of the Roche 
lobe ljVerbuntlll99.lt) for c ompanion masses of 4, 8, and 
16 M , we find C following ITben fc Tutukovlll996|) Vns ~ 
40, 560, and 1000 km s , with a possible small contri- 
bution from the prior systemic velocity (< 100 km s _1 ). 
Thus the disruption of a close binary is very unlikely to 
solely account for the birth velocity of B1508+55, and a 
natal (or post-natal) kick is indicated. 

Natal kicks from asymmetric S N explosions are 
believed to be ubiquitous (e.g.. iFrver fe Kalogera 
1997 1 iPortegies Zwart fc van den Heuvellll999HLai et al. 
2001), but the driving mechanism remains unclear. 
Global asymmetric perturbations during SN core col- 
lapse may produce h ydrodynamically driven kicks 
ijBurrows fe Hay esll996|). In recent two-dimensional sim- 
ulations l)Scheck et al J 120041) . velocities comparable to 
that measured for B1508+55 have been obtained. How- 
ever, the first full three-d imensional simulations of SN 
core collapse l)Frverll2004ll suggest that as the proto-NS 
moves through the convective region, asymmetric down- 
flows and neutrino emission damp out the initial kick 
velocities, producing final kicks < 200 km s _1 . Such 
kicks are insufficient to account for the observed velocity 
of B1508+55, but the issue is not yet settled. 

If convective instabilities in core collapse simulations 
fail to reproduce large kick velocities, this may suggest 
the existence of neutrino driven kicks, and require the 
presence of ex treme magnetic fields (B > 10 16 G) durin g 
core collapse ll.Tanka fc RaffeltJll999t lArras fc Lal ll999'l. 
Such a high B may also allow a post-nat al rocket effect 
to accelerate the pulsar (jLai et alJ [20011. However, a 
kick driven by the NS £?-field should produce alignment 
between the spin axis and the velocity vector due to rota- 
tion averaging. The spin axis may be inferred from the 
polarization sweep of the radio pulse for some pulsars 
(Rad hakrishnan fc Gookeill969l however, this is model 
dependent). For B1508+55, the pulse centroid polariza- 
tion angle (after correctin g for Faraday rotation a long the 
line of sight) is 161° ± 5° (peshpande et al.lll999|) . which 
differs from the proper motion angle by ~ 71° (or ~ 19° 
for emission in the orthogonal mode). Assuming it sur- 



vives future scrutiny, the measurement does not support 
alignment of the spin axis and Uns- 

We conclude that a kick is required to provide 
B 1508+55 with its high birth velocity, but that a single 
model is not uniquely preferred based on our current un- 
derstanding. It is possible that binary disruption and a 
natal kick act in concert for this object. Such a composite 
scenario would produce a misalignment between the spin 
axis and velocity vector, as is (possibly) observed, with 
the relative weights of the contributions depending on 
whether the emission is in regular or orthogonal modes. 
However, such a scenario is not well supported unless the 
polarization angle is verified and a run away binary com- 
panio n is identified in the future (e.g. iVlemmings et all 

2004) . Finally, we note th at in their seminal paper, 
iDuncan fc Thompson (1992) argued that magnetars had 
high space velocities ~ 1000 km s _1 due to kicks that 
required intense magnetic fields (B > 10 15 G), and that 
such strong kicks would not be accessible to ordinary 
pulsars (like B1508+55, with B = 2 x 10 12 G). If this ar- 
gument is to hold, then B1508+55 had a large birth mag- 
netic field that rapidly decayed. It is more likely that B- 
driven kicks are not required for high velocities, and that 
large natal kicks are not a defining difference between 
pulsars and magnetars. Rather, the dividing line between 
them may lie elsewhere in their birth a nd evolutionary 
histo ry (e.g., in their progenitor mass; iGaensler et alJ 

2005) ." 

Ongoing deep pulsar searches, and especially searches 
at high Galactic latitudes, coupled with long-term VLBA 
astrometry programs, are likely to discover further high- 
velocity pulsars and continue to elucidate the velocity 
distribution and birth circumstances of neutron stars. 
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